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Abstract 

We derive mass formulae for the ground state, = 0“ and 1“, and first excited even-parity, 
= O'*" and 1"*“, charmed mesons including one loop chiral corrections and 0(l/mc) counterterms 
in heavy hadron chiral perturbation theory. We show that including these counterterms is critical 
for fitting the current data. We find that certain parameter relations in the parity doubling model 
are not renormalized at one loop, providing a natural explanation for the observed equality of the 
hyperfine splittings of ground state and excited doublets. 
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I. INTRODUCTION 


Excited charmed mesons with angular momentum and parity = O’*" and I"*" have been 
observed in several experiments. The masses of the = 0+ and charmed strange mesons, 
Us(2317) and Ds{24:60) |l[ Q|, are below threshold for decays into ground state charmed 
mesons and kaons. The only strong decay modes are via isospin-violating 7r° emission, 


making the states quite narrow (T < 5.5 MeV). Other experiments claim to observe 

the nonstrange = O’*" and I’*' states. These states can decay to the ground states by 
S'-wave pion emission and therefore are quite broad (T ~ 300 MeV). 

The spectrum of the = 0^ and charmed mesons presents a number of puzzles 
for theory. Before their discovery, quark model and lattice calculations predicted that the 
masses of the = O’*" and charmed strange mesons would be signihcantly higher than 


observed 


Q 0, y, y, Q- 


Further, the hyperhne splittings of all ground state charmed 
mesons and the hyperhne splitting of the Ds{2317) and ^^(2460) are all equal to within 
2%. This is surprising because there is no obvious symmetry of quantum chromodynamics 
(QCD) which predicts these equalities. Finally, the SU{3) splittings of the = O'*" and I’*' 
charmed mesons are much smaller than theoretical expectations. 

In the heavy quark limit, the coupling of the heavy quark spin to the light degrees of 
freedom in the heavy meson vanishes and the angular momentum and parity of the light 
degrees of freedom, can be used to classify heavy meson states. The spectrum consists of 
degenerate heavy meson doublets with dehnite j^. The = 0“ and 1“ heavy mesons are 
members of the | ground state doublet. The lowest lying excited states, the = 0+ 
and heavy mesons, are members of the doublet. There is also an excited doublet 

of heavy mesons with = I"*", whose members have and 2~^. The = I"*" mesons 

decay to the ground state by D-wave pion emission, typically have widths T ~ 20 MeV, and 
therefore have well-measured masses. The hyperhne splittings for all of these heavy quark 
doublets are suppressed by l/mq, where mg is the heavy quark mass. 

The experimental data on the masses of the known charmed mesons is summarized in 
Table m The lowest lying havor SU{3) anti-triplets are = 0“ {cu, cd, cs) ~ 

Df) and = 1~ {D*^, The hrst excited states are = O’*" (Dq, DQ,DQg) 

and = 1+ {Df, Df', Df^). The members of the doublet are = 1+ (U°, 

Di^ D^g) and = 2^ (Ug, D^g)- Not shown is a narrow charmed strange meson. 
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jp JP 

cu 


cd 


cs 


3/2+ 2+ 

name 

M(MeV) 

2458.9 ± 2.0 

name 

Dt 

M(MeV) 

2459 ± 4 

name 

Dt2 

M(MeV) 

2572.4 ±1.5 

3/2+ 1+ 

D? 

2422.2 ± 1.8 

Dt 

2427 ± 5 

Dti 

2535.4 ±0.6 

1/2+ 1+ 

D?' 

2438 ± 31 

- 

- 

Dt{ 

2459.3 ± 1.3 

1/2+ 0+ 

nO 

-^0 

2308 ± 36 

- 

- 

Dto 

2317.4 ±0.9 

1/2" f- 

D*o 

2006.7 ±0.5 

D*+ 

2010.0 ±0.5 


2112.1 ±0.7 

1/2- 0- 

DO 

1864.6 ±0.5 

D+ 

1869.4 ±0.5 

Dt 

1968.3 ± 0.5 


TABLE I: The spectrum of charmed mesons, is the angular momentum and parity of the light 
degrees of freedom. is the angular momentum and parity of the meson. 


D^(2632), recently observed by the SELEX collaboration |ll|. The spin and parity of this 
meson and its place in the charmed meson spectrum is currently unknown. For all mesons 


except the nonstrange doublet, we use numbers from the Particle Data Group 


For nonstrange = 


1 + 


mesons, we use the Belle 


average the CLEO 


3 and Belle Q 


Q 


id- 


measurement of the Dq mass and 


measurements of the mass. 


As stated earlier, the hyperhne splittings of the ^ and doublets are nearly 


equal. The known hyperhne splittings of the ^ and charmed mesons are: 


1 + 


Tn£)*o — TTljyO 

mD*+ — mo+ 
mDi+ - m^+ 

-^sl -^sO 

mDOi — 


142.1 ± 0.07 MeV 
140.6 ± 0.1 MeV 

143.8 ±0.4 MeV 

141.9 ± 1.6 MeV 
130 ± 48 MeV. 


( 1 ) 


Here, t 
Group 


le hrst three numbers are the hyperhne splittings quoted by the 
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The last two numbers are obtained by taking the diherence 


Particle Data 
of the masses 


^ The FOCUS collaboration reports structures in excess of bacteround in the and invariant 

mass spectra which could be interpreted as scalar resonances However, if these resonances exist their 
masses are 99 MeV higher than the Belle measurement and 80 MeV higher than the mass of the It 
is implausible that such resonances are related to the by SU (3) symmetry so we do not use this data 
to determine the = 0+ nonstrange meson masses. 
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in Table lU The error in the last two lines of Eq. o is obtained by adding the errors in 
the individual masses in quadrature. All four hyperhne splittings which have been mea¬ 
sured accurately are ~ 142 MeV to within 2 MeV or less. Hyperhne splittings in different 
heavy quark doublets are unrelated by heavy quark symmetry. For example, the hyperhne 
splitting for doublets is ~ 40 MeV, which dihers signihcantly from the | 

and hyperhne splittings. In the SU{3) limit, the hyperhne splittings of nonstrange and 
strange ground state mesons are the same. That this SU{3) prediction holds to within 2% 
is surprising given the typical size of SU{3) breaking ehects in QCD. 

Another puzzling feature of the spectrum is the pattern of SU{3) violation in the splittings 
within the even-parity doublets. Finite light quark (m„, rrid, and m^) masses and electro¬ 
magnetic ehects cause havor-splitting among the mesons. The isospin splitting seen in the 
charmed meson mass spectrum is of expected size, but the splitting between the strange 
and nonstrange sector is unexpected. The mass diherence between strange and nonstrange 
mesons whose other quantum numbers are identical is expected to be ~ 100 MeV. For the 
ground state charmed mesons this is indeed the case. For the excited states, however, the 
SU(3) breaking is 


rriTj+i — mjyoi = 21 ±31 MeV 

^sl 1 

m^+ — m^o = 9 ± 36 MeV. (2) 


Even allowing for the large errors due to the uncertainty in the masses of the nonstrange 
jP = charmed mesons, the SU{3) splitting is far below theoretical expectations. 

The Ds{2317) and Zls(2460) are only 40 MeV below the DK and D*K threshold, re¬ 
spectively. This fact as well as the puzzles mentioned above have led to the hypothesis that 

Q,Q,Q- 


they are bound states of and K 


Several papers analyze the spectroscopy 


of excited charm mesons by extending the quark model to include couplings to the DK 


continuum. This couplec 
chiral quark models (l7 


channel effect has been analyzed within the quark model 


18| | as well as unitarized meson models jlOl. l20l. |2l|. The unita- 


Q. 


rized meson model has also been used to make predictions for the spectroscopy of excited 
B mesons 1^. However, if one assumes that the Ds{2317) and ^^(2460) are nonrela- 
tivistic DK and D*K bound states, respectively, heavy-hadron chiral perturbation theory 


(HHyPT) 


can be used to predict their electromagnetic branching ratios. These pre 


dictions are found to be in serious disagreement with experiment 


24| . On the other hand. 
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if one assumes that the Ds{2317) and ^*(2460) are conventional states, then HH^PT pre¬ 
dictions for strong and electromagnetic decays are consistent with available data 


An alternative interpretation of these particles as exotic csqq tetraquarks has also been pro¬ 


posed 
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For a review of theoretical work on ^^(2317) and 77^(2460), 


see Ref. 

In this paper we analyze the spectroscopy of charmed mesons using HH^PT. This theory 
can be used to analyze the low energy strong interactions of heavy mesons in a systematic 
expansion in light quark masses, rUg, and inverse heavy quark masses, XjrtiQ. Nonanalytic 


corrections from loops with Goldstone bosons can be calculated in this formalism, 
masses of the ground state heavy mesons have been studied in the heavy quark limit 


'he 


32 


33|, including leading corrections from hnite heavy quark masses and nonzero light quark 



37|, 


38, 


39|. These papers use a version of HHyPT which includes only 


masses |3J, |35|, |3 


the lowest lying ^ heavy quark doublets. Many recent studies of excited = O’*" 
and heavy mesons use Lagrangians that include only = 4 and = I"*" heavy quark 
doublets as explicit degrees of freedom. However, the excited doublets are only 

separated from the = I"*" doublets by ^ 130 MeV. Further, the doublets couple 

to the = 4’’“ doublets at leading order in the chiral expansion, while the coupling of the 
jp _ 1 + (^loublets to the ground state doublets is higher order in the chiral expansion Q. 

o I 

For these reasons, loops with virtual excited | could have important effects on the 
physics of = 4^ doublets. In this paper we will study the version of HHy;PT containing 
only the = 4 and heavy quark doublets and leave investigation of loop effects 

from more highly excited states for future work. 

A model of heavy mesons closely related to HHy;PT is the parity doubling model of 


Refs. 


E, 3, 


44| . The parity doubling model is the analog of the linear sigma model for 


heavy mesons. Heavy meson doublets transforming linearly under SU{3)l x SU{3)r couple 
in a chirally invariant way to a held S transforming in the (3,3) of SU{3)l x SU{3)r. The 
held S develops a vacuum expectation value and the resulting theory of heavy mesons has 
the same form as HHy;PT for the low lying odd- and even-parity doublets. Unlike HHy;PT, 
the parity doubling model predicts relationships among otherwise independent parameters 
in the theory. One important prediction is that the hyperhne splittings of the = 4 
and = I’*' doublets are equal at tree level. This interesting prediction could partially 
explain the observed pattern of heavy meson hyperhne splittings, but it is not clear from 
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Refs 


• H 


m 


43 


44| whether this prediction survives beyond tree level. This is a concern 


because loop corrections in HHyPT can be signihcant. 

In this paper, we calculate the one loop HHyPT corrections to the masses of | 
and jP = heavy meson doublets. We include all 0{l/mQ) heavy quark spin symmetry 
violating operators that appear to this order. A brief review of the HHyPT formalism is 
given in section II and explicit formulae for the masses at one loop appear in the Appendix. 
In section III, we attempt to £t the observed mass spectrum with our one-loop formulae. 
The large number of free parameters makes it possible to reproduce the spectrum of | 
and 7 ^ = charmed mesons. In the mn —>■ oo limit our calculation of the difference of 

n 

the SU( 3) splittings in HHyPT agrees with Ref. |4^. Our analysis differs from that in 
Ref. in that we include I/tuq operators and perform a global fit to the spectrum with 
all counterterms treated as free parameters. In the approximation used in Ref. there is 
a single counterterm constrained using lattice data. 

In section IV, we examine corrections to the hyperhne splittings and discuss the natural¬ 
ness of the parity doubling model. The parity doubling model predicts that the hyperhne 
splittings and the magnitudes of the axial couplings of the = 
are equal at tree level. We hnd that these parameter relations are preserved by the one loop 
corrections so that the model provides a natural explanation for the equality of hyperhne 
splittings. Finally, in section V, we use heavy quark ehective theory (HQET) to estimate 
the masses of the R mesons, which have not yet been observed. These predictions 

may be helpful to experimentalists looking for these states. 


and = \^ doublets 


II. HHxPT MASS COUNTERTERMS 

In HKyPT, the ground state doublet is represented by the helds 


= 


- Hals) , 


( 3 ) 


where a is an SU (3) index. In the charm sector, Ha consists of the (R°, R+, Df) pseudoscalar 
mesons and Rh are the (D*^, D*~^, vector mesons. The doublet is represented 


by the helds jT 


Sa = - Sa) , 


( 4 ) 
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where the scalar states in the charm sector are Sa = D^a and the axial vectors are = D'^^. 
The kinetic terms of these helds are included in: 

Jl^kinetic ^ _Tr[F,(zU ■ D,, - 6h Sab)H,] + Trp,(*U ■ D,, - 6s , (5) 

where 6h and 6s are the residual masses of the H and S helds, respectively, and is 
the chirally covariant derivative. In the theory with only H helds one is free to set 6h = 
0. Since the only dimensionful parameters entering the loops in this theory are hyperhne 
splittings and meson masses, the UV divergences (in dimensional regularization) vanish 
in the rriq —>■ 0 and tuq cxo limit. Divergences in loop corrections are canceled by 
counterterms which are 0{mq) or 0{l/mQ). Once the S helds are added to the theory, there 
is another dimensionful quantity, 6s — 6 h, which does not vanish as rriq —>■ 0 and rriq — »• oo. 
The H self-energy diagrams with virtual S helds give a UV divergent contribution which 
survives in the rriq 0 and mg —> cx) limit. Such a divergence must be canceled by a mass 
counterterm which respects SU{3) and heavy-quark spin symmetry and the only available 
counterterm is 6h'^^H alia- However, after one-loop divergences are canceled one is free to 
dehne the hnite part of 6h for convenience. 

The helds have axial couplings to the pseudo-Goldstone bosons, 

^ g + g' Tr+ h Tr+ h.c.] , (6) 

where g', and h are dimensionless constants to be determined from experiment. The other 
terms in the Lagrangian required are higher order mass counterterms. We use the notation 
of Ref. and generalize it to include the S helds as well as the H helds. 

+aH"^AHaHh] ml^ - asTr)^^^;,] -F aHTr[hf„if„] ml^ - cr5Tr[S'aS'„] ml^ 

ml + ml 

ml + ml , ( 7 ) 

o o 

where = l{^rnq^+^^mq^^)i,a, rriq = diag(m„, rrid, and ^ = a/S = exp(in//), where H 
is the matrix of Goldstone bosons. The hrst line in Eq. o contains spin-symmetry violating 
operators which give rise to hyperhne splittings. The second line contains terms which are 
0{mq) and preserve heavy-quark spin symmetry. Finally, the last two lines contain terms 
which are 0{mq) and violate heavy-quark spin symmetry. 
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FIG. 1: One-loop self energy diagrams for the H and S fields. H fields are single lines, S fields are 
double lines and Goldstone bosons are dashed lines. 


At tree level the residual masses are 


m 


Ha 


= 6h - + aHm + anma 


- A m 
4 


m 


1. _ Ia(o-) _ 1a 

= Oh + + <Jh m + ttH rua + ^ + '^^h 

= 6s- + asm + asma- ^A^f^ fn - ^A^“^ 


77 ? 


(a) 


H* 


mi 


mi. 


. 1a — 1a 

= ds + - As + asm + asma + - A^ 




m 


1 ^ 

+ ~Ag 

4 * 


(a) 


m„ 


( 8 ) 


where ma = {mu,md,ms) and m = mu + md + mg. Here the asterisk denotes the spin-1 
member of the heavy quark doublet. In the isospin limit mu = md- HHyPT is a double 
expansion in Aqcd/'^q and Q/A^, where Q ~ mT,,mK,fnri and A^ = Anf 1.5 GeV. The 
parameters 6h, 6s, Ah, and As are treated as 0{Q) in the power connting [2^. Since 
mq oc ~ the remaining terms in Eq. (jH} are formally higher order in the power 
connting. The loop corrections to the masses are shown in Fig. ^ Single lines represent 
the H helds and double lines represent the S helds. All diagrams are 0{Q^). The loop 
corrections are regulated using dimensional regularization. Complete one loop expressions 
for the masses are given in the Appendix. 


III. CHARMED MESON SPECTRUM 

In this section we analyze the charmed meson spectrum using the one-loop mass formnlae 
given in the Appendix. We will work in the isospin limit, where the masses of Hi and H 2 , 
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for instance, are identical. Then there are eight different residnal masses: rriHii 
'm-H*, ^Si, fnss, fnsi, and ms*. To determine the experimental valnes of and mn*, 
we average the masses of the two known isospin states. The residnal masses are defined 
to be the difference between the real masses and an arbitrarily chosen reference mass of 
0{mQ). We will measnre all masses relative to the nonstrange spin averaged H mass, so 
+ 3mH*)/4: = 0. Therefore, the experimentally measnred residnal masses we will fit to 

are: 


niHi = —106.1 MeV = —4.75MeV mn* = 35.4 MeV mn* = 139.1 MeV 

= 335.0MeV = 344.4MeV = 465.0MeV ms* = 486.3 MeV. (9) 

The tree level expressions in Eq. (jS)) reprodnce these valnes with 5s + (Tsm — 5 h — c^nm = 
432 ±26 MeV, = 146 ±1 MeV, A^VA^^^m = 129 ±50 MeV, an = 1.14±0.06, 

as = 0.21 ± 0.29, = -0.03 ± 0.01, and A^^^ = 0.14 ± 0.55. The errors nsed to 

obtain this fit are the experimental ones, dominated by the nncertainty in the nonstrange 
O’*" and 1^ masses. This gives rise to the large nncertainties seen in parameters in the 
Lagrangian involving the S helds. The hts presented in this section nse Mathematica Q 
and/or Miniiir 

The loop corrections depend on eleven parameters: g, g', h, an, cis, A^\ 5H + crHm, 

5s + '^sm, Ah + m, and As + A^g'^rn. The parameters an, 0's, A^\ and A^^^ cannot be 
separately determined becanse they always appear in linear combination with the parameters 
5h, 5s, Ah, and As, respectively. Below we will absorb the contribntion of the parameters 
an, 0's, Ah\ A^^ into the measnred valnes of 5h, 5s, Ah, and As, respectively. 

An analysis of D* decays using a one-loop calculation without explicit excited states 
yields g = 0.271°;°® j^. From the widths of the nonstrange resonances observed by Belle 


we have extracted h = 0.69 ± 0.09 at tree level 


24| . Both couplings are of order unity and 


therefore consistent with naive power counting. The remain ing parameters are unknown. 

We use / = 120 MeV, which is the value extracted in Ref. ^ using the one loop formulae 
for pion and kaon decay constants, first derived in Ref. [^. We set = md = 4 MeV 
and mg = 90 MeV. Below we show several different fits. In the first case we fix g and h to 
the values (given above) extracted from previous analyses. This leaves nine remaining free 
parameters. Performing a least chi-squared fit to the meson spectrum, using experimental 
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uncertainties, we obtain the following central values 


tuhi — — 106MeV = —5MeV rnn* = 35MeV Triu* = 139 MeV 

= 160 MeV ms^ = 344 MeV ms* = 296 MeV ms* = 486 MeV. (10) 

The parameters extracted from this fit are: g' = 0.09± 0.03, 5h = —83±3 MeV, 5s = 244±1 
MeV, Ah = 133 ± 2 MeV, A 5 = 136 ± 1 MeV, an = 1-70 ± 0.01, as = 0.25 ± 0.08, 
Ah^ = —0.07 ± 0.01, and A^^ = 0.04 ± 0.03. Six of the mass parameters are reproduced 
quite well while ms^ and ms* are lower than the central values of experiments by about 
175 and 169 MeV, respectively. This qualitative picture persists without much sensitivity 
to the value of g'. However, these fits are not very good and such a procedure may not be 
very realistic. The values of g and h used above were extracted using a £t to a one-loop 
calculation not including the S fields, and a tree-level fit, respectively. There is no reason to 
believe that these values are the ones which are appropriate for a calculation that includes 
graphs with internal S states. Note that large changes between tree- versus loop-extracted 
parameter values do not necessarily indicate poor convergence; what is important is that 
the observables do not suffer large changes between orders. 

If we include g and h as free parameters in an 11-parameter fit, there are many solutions 
which yield central values identical to the experimental residual masses given in Eq. Q. In 
addition to the experimental errors we also include 20 % “theoretical” errors to mimic the 
fact that our analysis is only accurate to 0{Q^). The masses obtained are then accompanied 
by errors at the ± 30 to 40 MeV level. Examples of parameter sets which give these results 
are: 

(a) \g\ = 1.15±0.06, \g'\ = 0.90T0.06, \h\ = 2.3T0.2, 5h = 195±41 MeV, 5s = 332±31 

MeV, Ah = 465 ± 24 MeV, Ag = 597 ± 28 MeV, an = 7 ± 1, as = -4 ± 1, = 

-4.4 ± 0.7, and A^“^ = -10 ± 2. 

(b) 1^1 = 0.65±0.06, \g'\ = 0.89±0.08, \h\ = 0.2±0.1, 5h = 117±21 MeV, 5s = 646±40 
MeV, Ah = 68 ± 42 MeV, As = 447 ± 23 MeV, an = 3.8 ± 0.7, as = 3.1 ± 0.7, 
A^^ = -0.3 ± 1, and A^“^ = -2.8 ± 1. 

(c) 1^1 = 0.89T0.07, l^'l = 0.24T0.13, \h\ = 0.98T0.11, 5h = 203±39 MeV, 5s = 399±26 
MeV, Ah = 242 ± 25 MeV, As = 116 ± 59 MeV, an = 5.8 ± 1.1, as = -1.4 ± 1.5, 
A^h^ = -1.7 ± 0.8, and A^ = 2.1 ± 1.7. 
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(d) 1^1 = l^'l = 0.70 ± 0.03, \h\ = 2.4 ± 0.2, 5h = 114 ± 64 MeV, 5s = 231 ± 36 MeV, 
Ah = 682 ± 4 MeV, an = 4.3 ± 0.7, as = -3.0 ± 2.1, = -0.89 ± 0.96, and 

= —2.7 ± 0.9. (In this fit, the constraint A 5 = Ah + 30 MeV was used.) 


There are clearly many local minima which Minuit m may hnd. Some of these values of 
g and h signihcantly exceed values extracted from experiment in Refs, 
respectively. They also exceed estimates based on the quark model 


49 


lattice QCD simulations 


53 


as well as sum rule constraints 


3 - 


and Ref. 


m, 


m, extraction from 
Again, however, it 


is not clear what can be concluded when comparing parameters which are by themselves 
unphysical and whose dehnition depends upon the details of a calculation. Of perhaps more 
concern is that these hts produce large values for the hyperfine coefficients. The operators 
which cause hyperhne splitting should be I/tuq suppressed compared to the leading order 
ones. Set (d) is an example of a solution where is near l^f'l and As is within 30 MeV 
of Ah- The relevance of that result will become apparent in the next section. Finally, an 
example £t where g and h are restricted to lie between 0 and 1 yields central values of 
parameters as follows: 


\g\ = 0.43, = 0, h = 0.31, 6h = 25 MeV, 6s = 443 MeV, 

Ah = 124 MeV, A 5 = 131 MeV, an = 2.4, as = -0.3, A^^^ = -0.2, A^“^ = 0.1. 


These parameter values lead to a prediction for the mass spectrum that also agrees with 

Eq. ©. 

The under const rained nature of the various fits makes strong conclusions impossible. 
In particular, the uncertainty in the parameter space is very large and the uncertainty in 
individual parameters is much greater than indicated by the errors quoted in the individual 
hts listed above. The situation should improve with a global £t to both masses and decay 
rates which uses a consistent set of next-to-leading order calculations that include the excited 
states. This work is in progress. 


IV. HYPERFINE SPLITTINGS 

In this section we study the one loop corrections to the hyperhne splittings to see if 
HHy;PT can provide insight into the observed near equality of the hyperhne splittings. Using 


11 










the formulae in the Appendix we find that the next-to-leading order difference between even- 
parity and odd-parity hyperhne splittings in the strange sector is given by 


{ms* -nis^) - {mm - mn^) = {m% -m%)- {m^H: - 


g'^ 12 


P L3 
4 


^24 

+ gKi{m°s^ - m^) + -Ki{ml, - m^) 


P 


^-Ki{0,mr,) -2Ki{m%^ -mgi^^K) - -Ki{msi 
2 '"2 2 4 


9 

4 2 

+ gKi{0, mr,) -2Ki{m%. - m^H^^mK) - -Ki{m'lj, 

2 

2K2{mH* - m^sp'rriK) + -K2{m^H* “ ~ 

2 2 
- o^2(m^3 - ml^,mrj) - 2K2{m%, - mH*,mK) - -K2{mm - m^) 


3 

,0 _o 


P 


+ 2 K 2 { m^Si - Ph 3 , ^ k ) + 2^:2(77153 - 771 ^ 3 , m^) 


( 11 ) 


Suppose one imposes at tree level the condition that all hyperhne splittings in each of the 
doublets are degenerate: 


m 


H* 


m 


Ha 


= ml 


m 


Sa 


= A 


( 12 ) 


This can be arranged by invoking the tree level prediction of the parity doubling model, 
Ah = As = A and neglecting the terms proportional to mq in Eq. (jHl), which are formally 
higher order in the power counting. Then m^^, — 7?7^* = m^^ — m^^^ and it is easy to verify 
that all contributions proportional to vanish, and the remaining terms are: 


v 2 


( 7715 . - 77753 ) - {^H* - mHa) = 


P l3 


2 16 

A'i(-A, 777^) -h -A'i(-A, 777^) -—A^i (0, 777i^) 


9 


9 


2.1 12 


-2A'i(A, 777i^) - -iri(A,7?7^) 




16 2 
-hYA'l(0,777i^) - 2A'i(A,7?7i^) - -A'i(A,7?7^) 


(13) 


This vanishes if = g'^, which is consistent with the parity doubling model prediction. A 
similar cancellation occurs for the nonstrange hyperhne splittings. So the parity doubling 
model explanation for the equality of the = 1 and P hyperhne splittings is robust in 
the sense that one loop corrections do not spoil the prediction. 

The parity doubling model prediction for the axial couplings and hyperhne splittings sin¬ 
gles out a subspace of the parameter space of HH^PT that is preserved under renormalization 
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FIG. 2: One-loop diagrams for renormalization of the coupling g. H fields are single lines, S 
fields are double lines and Goldstone bosons are dotted lines. Diagrams for renormalization of the 
coupling g' are obtained by interchanging H and S fields. 


group evolution. From our mass formulae it is easy to derive the following renormalization 
group equations for the renormalized parameters Ah and A^: 


d . 


-t—^Ah = 


■A' 


97r2/2 ^ 


(14) 


h 2 


/i 


2 d o 

2 Ac = — —A^ 


:{As-Ah) 


3(55 — duY — — Ah){5s — 5h) + — AhY 

2 ib 


djj? 


97r2/2 ^ 

~*~ 37 ^ 2 p ~ 


(15) 


3(^5 — 5hY — 7^{As — AH){ds — Sh) + ~ AhY 

2 lb 


which leads to 


. ^(/a|-/A3,) 


(16) 


2 /2,2 

+ - Ah) 


3(5s — — -{As — Ah){Ss — Sh) + tw('^s' — Ah)"^ 

2 lb 


We also derive the one loop renormalization group equation for the couplings g and g'. 
For this we need the wavefunction renormalization of the helds H and S', which is obtained 
from the graphs in Fig. Q and the one loop corrections to the axial couplings. The relevant 
graphs for the renormalization of g are shown in Fig.|2l and the graphs for g' can be obtained 
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from those in Fig. |21by interchanging H and S lines. Note that we only need the ultraviolet 
divergences of these graphs to obtain the renormalization group equation. Furthermore, 
the counterterms for the wavefunction renormalization and the axial couplings are dehned 
to be independent of rriq and tuq. ^ Ultraviolet divergences proportional to rUg and I/tuq 
are absorbed into higher order counterterms. For example, a divergence proportional to niq 
in the one-loop correction to the axial coupling of the H helds should be renormalized by 
counterterms with structures like TT[HaHb^bc'j 5 \ml^, TrlHaHh^balbl^ic^ etc. Therefore we 
can ignore Goldstone boson masses and hyperhne splittings in computing the ultraviolet 
divergences, which greatly simplihes the calculation. The graphs in Figs. EK, andEfe 

vanish in this limit because the integrals are scaleless. Graphs in Figs. Et and Eli do not 
contribute either. This is because the H-S-n coupling in Figs.Et andEti gives a factor oiv-k, 
where is the four-momentum of the external Goldstone boson. Ultraviolet divergences in 
Figs. Efc andEti are proportional to v-k and are canceled by counterterms with an additional 
covariant derivative acting on the helds such as TY[HaHbiv ■ Dbc^-ca'Jb]- Therefore, all 
that is needed to obtain the running of g are the ultraviolet divergent parts of Fig. ^ and 
Eb in the limit where Goldstone bosons and hyperhne splittings are neglected. The running 
of g' is obtained from Fig. ^ and the analog Fig. Eb with S and H lines interchanged. The 
result can be obtained from the corresponding graphs for the renormalization of g by simply 
substituting g ^ g' and Ss — Sh ^ — ^h)- The renormalization group equations for g 

and g' are 

which can be rewritten as 

- g') = - g') • (18) 

To understand the signihcance of this result, consider the naive quark model prediction 
g' = g/^ 1^ . From the renormalization group equations in Eq. m one sees that g and g' 

^ If the theory is not renormalized this way, dependence on the underlying theory parameters mq and toq 
would no longer be explicit in the chiral Lagrangian. 
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vary with changes of the renormalization scale in snch a way that the condition g' = g/3 can 
only hold at one valne of g,. The qnark model prediction is meaningless beyond tree level 
withont also specifying a particnlar renormalization scheme and scale at which the relation 
is expected to hold. However, ii g = -^g' holds at any /i, it will hold for all g (at least at one 
loop order). Also, if g"^ = g'"^ and A 5 = the right hand side of Eq. (fTT)|l vanishes. Thns 
the predictions of the parity donbling model, Ah = A 5 and g = —g', are invariant nnder 
renormalization gronp flow in HHyPT to one loop order. 


V. HQET AND PREDICTIONS FOR EXCITED B MESONS 


In this section, we comment on the theoretical expectations for the spectrnm of excited 
even-parity bottom mesons which have yet to be discovered. Onr HHyPT resnlts for the 
charmed meson spectrnm may be nsed, bnt there are nnknown 0{l/mQ) effects which 
make it difficnlt to obtain precise predictions for the B meson. For finite qnark masses, to 
obtain the bottom meson spectrnm from the charmed meson resnlts, the hyperhne operators 
shonld be rescaled by rric/mf,, which is not very well determined. Other parameters can 
receive 0 (Aqcd/''^c — AgcD/^^fe) corrections. For instance, the rednced kinetic energy of 
the b qnark signihcantly rednces the mass splitting of the H and S donblets in the b sector 
relative to what is observed in the charmed system. These 0{l/mQ) corrections introdnce 
signihcant nncertainty in HHyPT predictions. 

We will instead nse the Od/mo) HQET formnlae for the mass of a heavy hadron X 


which contains a heavy qnark Q 

= mg -1- A" 


Af Af 

+ nj-^ 


2m, 


Q 


2m, 


Q 


( 19 ) 


where and Xf are hadronic matrix elements of the HQET operators higDYh and 
gsha^‘'Gfj,yh/2, respectively, and nj = -|-1 for J = 1 states and nj = —3 for J = 0 states. 
The hrst l/mg correction, — A^/(2mQ), is the kinetic energy of the heavy qnark. The 
second l/mg correction contribntes to hyperhne splittings. The difference between the 
spin averaged masses of the = A and = i"*" mesons, = (3m^2 -|- m^^)/4 and 

-f- respectively, is given by 

= A^ - A^ - ^ ^ , (20) 

^ ^ 2mQ 2mQ' ^ ’ 
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which leads to the following formulae for the splitting of the even- and odd-parity states in 
the bottom sector: 



A recent global £t to B decays yields Af = —0.20 ± 0.06 GeV^ |59|. The parameter Af is 
unknown. From the spectroscopy of excited = |^ iA and B mesons, Ref. 0 extracts 
Ai^^ — Af^ = —0.23 GeV^, where A^^^ is the Ai matrix element for the doublet. The 

sign here indicates that the kinetic energy of the heavy quark in the excited heavy meson 
is larger than that in the ground state, which agrees with intuition. We expect the kinetic 
energy of the heavy quark in the = 2 states to be comparable to that of | states. 
To estimate with conservative errors, we take Af — Af^ = —0.2 ± 0.1 GeV^, rric = 1.4 
GeV, and mb = 4.8 GeV to hnd 


m^^ — = m'g — — 50 ± 25 MeV (22) 

In the bottom nonstrange sector, m^| = 5279 MeV and m^l = 5325 MeV, which yields 
m^l = 5314 and therefore Eq. (1^ predicts m^^^ = 5696 ± 30 ± 25 MeV. The hrst error 
comes from the uncertainty in the charm nonstrange masses and the second is 

the estimated uncertainty in Af. These states are well above the threshold for S'-wave pion 
decays to the ground state and should be broad like their charm counterparts. 

In the bottom strange sector, only the 0“ state with mass m^^ = 5370 MeV has been 
observed. To proceed we need to estimate the mass of the bottom strange 1“ state. Note 
that 


m 


(b) 

H* 


m 


ib) 

H 


(b) 

moi — m 


(b) 


m 


(c) 


(c) 


m 


(c) 


m 


(c) 


me 

mb 


(23) 

— m'fj' m'gi — m'g' 

up to 0{l/mQ) corrections. Thus, all the hyperhne splittings in the bottom sector are related 
to those in the charm sector by a universal factor. Gombining this with the measured 
value of m^l — m^| leads to the prediction that m^l — m^^ = m^* — m^^ = 46 MeV, 
and m^gl — = 42 MeV. These predictions have approximately 25% uncertainty due 

to higher order 0{AQc-D/mc — AgcD/’^fe) corrections and the prediction for mgl — mg^ 
estimate has an additional 20% uncertainty due to the poorly known m^* and m^g^ masses. 
Given these hyperhne splittings, one expects m^H^ = 5404 MeV and then Eq. (j22j) predicts 
m^g = 5702 ± 25 MeV. Here the error is dominated by our ignorance of A^. Note that 
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the excited bottom strange mesons are expected to lie well below the threshold for decays 
to ground state B mesons and kaons and should be narrow like charmed strange 

mesons. 


VI. CONCLUSIONS 


We have enumerated the leading and subleading operators which describe the even-parity 
charmed meson masses in heavy hadron chiral perturbation theory (HHyPT). We performed 
a loop calculation to analyze the lowest lying even- and odd-parity charmed meson masses 
to 0{Q^). There are nominally eleven unknown parameters in the prediction, and only eight 
experimental masses. Two of the parameters, the axial coupling g for the lowest doublet 
of charmed mesons, and the coupling h which dominates the strong decay between the 
even-parity and ground state doublets, have been extracted from previous calculations. See 
Ref, Q and Ref, Q. ^ecfively, Howeven, fhe even-parity etatee were not tneluded tn the 
extraction of g in Ref. Also, the extraction of h was only performed at tree level. Since 
these values for g and h were not obtained under the same conditions as the mass calculations 
performed in this paper, it is not clear that the values should be used in our £t. Indeed, if 
the values from Refs. and Q are used, it is not possible to obtain the nonstrange even- 
parity masses as large as they are observed to be. If the g and h parameters are not hxed but 
simply constrained to he between 0 and 1, which is the prejudice from other analyses, then a 
£t to the even-parity masses is possible. Because of the numerous undetermined parameters, 
HHyPT can accomodate, but not explain, the unusual pattern of SU{3) breaking observed 
in the excited charmed meson spectrum. 

If we perform an unconstrained £t to the charmed meson mass spectrum using all eleven 
parameters, many solutions are possible, including ones whose values of g and h are not 
unreasonably far from their previously extracted values. However, then the parameter values 
found for the hyperhne operators are of concern. These hyperhne operators should have 
coefficients which scale as 0{AQ^^/mQ) relative to the 0{Aqcd) coefficients of the leading 
operators. The fact that unconstrained global fits find coefficients which are sometimes larger 
for the hyperhne operators than for the leading order operators may signal a problem in the 
l/rnq expansion. On the other hand, this may simply be a consequence of not properly 
incorporating the constraints on g and h from the decay widths. Before more definitive 
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statements can be made, a global fit including even-parity intermediate states and terms up 
to 0{Q^) for both the odd-parity and even-parity meson decay rates must be done. That 
will be the subject of a subsequent paper. 


Next we consider the parity doubling model introduced in Refs, 


42 


m 


While 


the parity doubling model is not a result of QCD, but requires additional assumptions, it 
is interesting because it provides an explanation of the observed equality of the hyperhne 
splitting in the even-parity doublet and the hyperhne splitting in the odd-parity doublet. 
QCD symmetries alone do not dictate any relationship between these hyperhne splittings. 
While the parity doubling model provides an explanation for the equality of the hyperhne 
splittings, the question we address here is whether it is a natural explanation. That is, does 
it survive beyond tree level? Is it stable under RG how? We hnd that there are “hxed 
lines” at = \g'\. (These are axial operator coefficients from Eq. (jU)).) That is, if at any 
time in their evolution g = g' or g = —g\ RG analysis shows that the relationship will be 
maintained. This in turn assures that if at tree level the parameters Ah and in Eq. o 
are equal, they remain so to one loop. This lends credence to the parity doubling model. 
The stability found in the parity doubling model does not exist for other models, such as the 
nonrelativistic quark model, which predicts g' = g/?). Going back to the parameter ht, we 
do hnd that solutions with \g\ near l^f'l are possible, as are hts with Ah near As- However, 
such hts yield values for Ah and As which are larger than expected by power counting. In 
addition, there are hts which reproduce the observed hyperhne splittings without ~ \g'\. 

Finally, we discuss how the charmed meson spectrum results can be used to make pre¬ 
dictions for the analog B meson spectrum. It is necessary to know the charm and bottom 
quark masses in order to rescale the operators, which brings in signihcant uncertainty. Also, 
there are additional l/mq operators with unknown parameters. However, it is possible to 
use heavy quark ehective theory to estimate that the even-parity strange spin-zero B meson 
has mass ~ 5667 MeV while its spin-one partner has mass ~ 5714 MeV. This places them 
below the threshold for decay to a kaon and the ground state B. Therefore, we expect 
narrow B* meson analogs to the narrow D* excited mesons. 

R.P.S. and T.M. are supported in part by DOE grant DE-FG02-96ER40945. T.M. is 
also supported in part by DOE grant DE-AG05-84ER40150. T.M. would like to thank the 
Aspen Genter for Physics and the KITP in Santa Barbara where parts of this work were 
completed. 
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VII. APPENDIX 


We express our results in terms of the functions 


Ai(r/,M) = 


K2iV:M) = 


( \ / 77 \ 

(—27]^ + 3M^?7) In —— ) + 2?7(?7^ — M^) F i — ) + Arj^ — brjM'^ 

\ /i^ / \MJ 

{—2rf + M’^rf) In [ ) + 2rf F (—^ + Ar]^ — t]M^ ( 

\ /i^ / \MJ 


where 


F{x) = 2- 


— x'^ \ti 


-Tan ^ 

2 




= —2 —- In (x + Va;^ — 1 

/v» \ 


bl < 1 


\x\ > 1 


The function Ki{r],M) appears whenever the virtual heavy meson inside the loop is in the 
same doublet as the external heavy meson, while -^’ 2 (^ 7 , M) appears when the virtual heavy 
meson is from the opposite parity doublet. 

In the limit M << 77 these functions become 


Aib,M) = 


K2{v,M) = 


-2,3 lnfl^U3,M^ Inf^ 


3M^ , /M 2 
+ In — 


4 , \ 4, 


-2, Inl^l+^M + 


4,2 \ 


In these equations we have dropped polynomials of ,, M. The functions Ki{r],M) and 
K 2 {r],M) have well-dehned M —^ 0 limits. Furthermore, the dependence on M is analytic 
when M/t] 0, so in this limit the S helds can be integrated out and their effect on the 
chiral corrections can be absorbed into local counterterms as expected. This limit is not 
relevant to the real world as , ~ M. In the opposite limit, , = 0, which is relevant for loops 
in which external and virtual heavy mesons are the same, 

KMM) = 

A 2 (,,M) = In +0(,3) (28) 

Including the one loop diagrams we hnd: 

^2 r-g g 

,2 pd g 

+ m,,) + F: 2 (m ^3 - • (29) 



mn, = niH^ + 


.0 , ^ 
/ 
r 


2 r 2 

- m%^, m^) 


P 


2K2{ms, - mH^,mK) + ^^ 2(^53 - 771 ^^, 771 ^ 


( 30 ) 


m/f* 


0 .9^ 

+ 753 L2 


2 1 r 3 I 

- m'jj*,77iP + -Ki{77i^jj^ - 


2 r3 1 

723 -i^'i(0,m^) + Ki{771h* - tti^.^ttik) 

r3 


P 


-K2{771s* - 771h*,77iP + -^^ 2 (^ 5 * - 771^^,771^ + K2{77ls* - TTli^) 


(31) 


0 , ^ 1 
+ —- 


2 1 r 9 

0 ^0 ^ I 7^0 ^0 


2Ki{771h^ - 771H*,771k) + -i^l(m^3 - 771 ^*,771^ 


— ^ 2Ki{ripj* -77i]j.,mK) + ^iri(0,m^) 

7^2 r 2 

-— \ 2K2{771%. - vPh* , ^x) + ^K2{771%. - 771%, , 771^) 


(32) 


ms:, = 


0,9 


/2 r 


3 1 

-Ki (m^* - 771% , TTlP + -^Ki ( 772 % - 771% ,77lP +Ki (771% - 771% , TUk) 


/^2 1-3 

-— y-K2(77l%^ - 77l%^,77lP + -K2(m%^ - 771%, 771% + K2(77l%^ - m^^,mx) 


.(33) 


mss = ms, + ^2 


,0 , ^ 

/ 


a 


2Kx(m%, - 771%^, 771k) + -Ki(77ps, - m%,,m% 


P 


2X2(771^, - ms,,mK) + -K 2 ( 771 k, - ms,,m% 


(34) 


ms* 


0 . 9 '" I 

+ 73 
.7^2 


-Ki{m% - m%,m% + -Ki{m% - m%,m% + Ki{m% - m%,mK) 


3 1 

P 3 ^ -Ki(0,m% + Ki(m% -m%,mK) 


P 


-K 2 (m%, - m%,m% + -K 2 (m%, - m%,m% + K 2 (m%, - m%,mK) 


.(35) 


ms* = 


0 , 9 '" I 

^s* + J,-, 


2Ki{m% - m%,mK) + -Ki{m% - m%,m% 


g%%2_ 

ps 
'P 


2Ki{m% -m%,mK) + -Ki{0,m% 

1 <5 


2K2(mH* - ms*,mK) + -i^2(m^3* - ms*,m% 


(36) 
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We agree with Ref. Q for the H helds in the limit where > 0, and 

rj/M -C 1. Our answer also agrees with that of Ref. which computes mass corrections 
to the H and S masses including SU (3) breaking corrections but not hyperhne splittings. 
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